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We report a simplification in the large N matrix mechanics of light-cone matrix field theories. 
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which exhibits an infinite number of new hidden conserved quantities at large N. 
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1 Introduction 



Light-cone quantum field theories [1-5] have been studied since 1949 and large N methods [6, 7] since 
1974. See Refs.[8, 9] for more complete referencing of the two developments. The idea of combining 
light-cone quantization with large N was pursued in [6, 10] in the context of summing planar 
Feynman diagrams, and in [11, 12] in the context of solving the large N energy eigenvalue problem 
in the color singlet sector of the light-cone Fock space. Here we study a different combination of 
the two developments-which leads to an algebraic reformulation of the light-cone theory at large 
N. 

In this paper we will focus on large N phase space methods [13-16] and in particular on the 
development known as large N matrix mechanics [16-18], [9], [19, 20]. In large N matrix mechanics, 
the Hilbert space is spanned by the singlet ground state and the dominant adjoint eigenstates which 
saturate the traces of the theory at large N, and the output is the so-called reduced formulation for 
both large N action [19] and Hamiltonian [9, 19] theories in terms of free algebras. The simplest 
free algebra found at large N is the Cuntz algebra, whose history in mathematics and physics can 
be found in Ref. [9]. In fact, a broad variety of free algebras has been observed [16], [9, 19] in large 
N matrix mechanics, including the Cuntz algebra, the symmetric Cuntz algebra, the interacting 
symmetric Cuntz algebras and infinite dimensional free algebras, as well as fermionic versions and 
Cuntz superalgebras. Moreover, the algebraic formulation has led to the discovery of new hidden 
conserved quantities [9] which appear only at large N, as well as various free-algebraic forms of the 
master field (see also Ref. [21]), and a numerical approach [20] to computations in many- matrix 
models at large N. 

But there has always been an enigma, called the opacity phenomenon in Ref. [9], in the large 
N matrix mechanics of equal-time quantized Hamiltonian formulations: Except for oscillators [9] 
and one matrix models [17], the reduced free-algebraic Hamiltonian of a given large N system is 
difficult to obtain in closed form. 

We report here an important simplification in the large N matrix mechanics of discretized 
light-cone matrix quantum field theories, including light cone models which are both local and 
non-local in x~ . The non-local models, including for example the number-conserving truncations 
of local models, can be considered as quantum-mechanical models in their own right. Because of 
the well-known absence of pure annihilation and pure creation terms in any light-cone Hamiltonian 
H. = P.~ we are able to obtain the explicit form of the large N reduced Hamiltonian entirely in 
terms of Cuntz operators. 

This opens up a free-algebraic playground of new reduced models, all of which exhibit an unex- 
pectedly simple algebraic structure. In particular, each of these models has new hidden conserved 
quantities at large N and the free-algebraic eigenvalue problems of the reduced Hamiltonians are 
surprisingly simple. 

The basic tool which we develop for the study of these models is the infinite dimensional 
algebra of all normal-ordered Cuntz operators, and this algebra also leads us to a special subset of 
the number-conserving models which apparently becomes integrable at large N. In particular, each 
of these models exhibits an infinite number of new hidden conserved quantities at large N, and 
moreover the eigenvalue problems of these models exhibit oscillator-like spectra. 

2 Light-cone Matrix Quantum Field Theories 

We begin with some general remarks about light-cone quantization of matrix field theories, deferring 
consideration of large N matrix mechanics to Section 3. 
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2.1 The Light-Cone Field 

The matrix field theory of most obvious physical interest is non-abelian gauge theory. But to simply 
illustrate the main ideas, we shall focus here on a much simpler example, the hermitian scalar N x TV 
matrix field (j) rs in two space-time dimensions. We define light-front coordinates = (x°±x 3 )/\/2, 
choosing x + as time, and referring the x~ coordinate to its conjugate momentum labelled by p + > 0. 
The field 4> rs then has the expansion 



4> rs (x ,X + ) 



(a(p + ,x + ) r 



-IX p n 



+ al(p + ,x + ) rs e A 



(1) 



with r, s = 1, . . . , N. The quantum nature of the field is fixed by imposing the equal time (equal 
x + ) commutation relations 

a(p + ) rSl a t (q + )tu 



SstS ru S(p + ~ q 



[a(p + ) rs ,a(q A 



)tu 



a){p + ) rs ,a)(q + ) tu 



0. 



(2) 



In the following we shall suppress the x + argument in expressions at equal time. A simple conse- 
quence of these commutation relations is the well-known equal time light-cone field commutation 
rule 



[4>vs{x ),4>tu{y )] 



-^ st 5 ru e(y - x ) 



(3) 



where e(x) = x/\x\, the sign of x. In particular, the fields at different spatial points do not commute, 
even at equal time. 

In our discussion we find it convenient to deal with discretely labelled operators so we discretize 
p+ = mS, m = 1,2, Putting a(mS) — ► a m /y/S, the commutation relations then read 



(fl"rn)rsi ifln)tu 
[("m)rsi {Q"n)tu\ 

(a m ) rs |0.) 



&mn 3st&ru 



h m)rs 



> (4) 



/ ti 



(.0|(a 



mJrs 







where |0.) is the Fock vacuum, and the expansion of the field reads 



<t>rs{x ) 



OO 



— ix m5 



+ (aln)r S e 



-ix mS^j 



-ix m5 



+ (fll)r S e 



-ix mS^j 



(4) 
(5) 

(6) 
(7) 



where K is a large integer serving as a cutoff on the high p + modes. With m an integer, we should 
restrict — ir < x~5 < n. The expansion (7) can be inverted 



iP"rn)rs ~ 

We also need the contraction factor 

dp + 
47rp+ 



tt/8 



7T J-tt/S 



D +ix mS 



4>rs{x ). 



I 



m=l 



K 



m=l 



Airm 



(8) 



(9) 



which arises in the process of normal ordering. 
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2.2 Partial Normal Ordering 



For the purposes of this paper the most important feature of light-cone coordinates is that the 
creation operators a){p + ) always create particles of positive p + . As seen in Eq. (1), we will assume 
that p + = is excluded, either by a small low momentum cutoff e, or in the context of discretized 
p+ = md, by the exclusion of m = 0. Only annihilation operators carry negative p + . Therefore the 
only p + conserving terms that can appear in the Hamiltonian contain both creation and annihilation 
operators, i.e. there can be no terms Tr a n or Tr a^ n , n > 0. As we shall see, this feature is 
responsible for a dramatic simplification in the free-algebraic description of the large N limit. 

Since the creation and annihilation operators are N x N matrices, as well as operators in 
Hilbert space, the process of normal ordering each term in the Hamiltonian can lead to awkward 
color descriptions. For example, the normal-ordered form of Tra^aa^a would be 

Tr : a)aa)a := a\ s a\ u a s ta ur . (10) 

Here there is a clash between the ordering for matrix multiplication and the ordering for Hilbert 
space operator multiplication on the r.h.s. Thus we refrain from completely normal ordering such 
terms. 

But we can partially normal order the operators in a trace, so there is always at least one 
annihilation operator on the extreme right or at least one creation operator on the extreme left 
without disturbing the coincidence of matrix ordering with quantum mechanical operator ordering. 
For example, a partially normal-ordered form of Trai a\a^a\ would be Tra\aia\a^. Such reorderings 
preserve the cyclic ordering of the operators within each trace. Note that in general a partially 
normal-ordered form is not unique, for example Tra^aia^a^ ^ Tio\<iz a \ a i although the operators 
are each in the same cyclic order. We shall see in Sees. 4 and 5, however, that the non-uniqueness 
disappears for N — > oo. Of course, in field theory the Hamiltonian is typically presented in a 
(local) form which is not even partially normal-ordered. But by exploiting the known commutation 
relations, we can always rearrange the operators in the partially normal-ordered form just described. 
Note that if there is at least one creation and one annihilation operator in a trace, partial normal 
ordering will lead to a generic structure 

Y ^{aln A m n a n ) (11) 

m,n 

where A is any matrix function of a p , aj, not necessarily itself normal-ordered. 
As a concrete example we consider the [/(iV)-invariant light-cone Hamiltonian 

H. = Pr = ^J dx'Txcj? + ±j dx-Ti^ (12) 

of a local two dimensional scalar field theory. Then the commutation relations (4) are used to 
obtain the partially normal-ordered form 

K K 
H. = E + r R(4 n A mn a n + ]T Tr(a m )B mn Tr{a n ). (13) 

m,n=l m,n=l 

The explicit forms of the quantities appearing here are 

E = [L 2 (iV 2 + l)A/4 + iV 2 L/2]^, B mn = 5 mn ^^- § (14) 
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A mn — C mn + D mnm2ms o J \ Jl2 cb\ n ^ + E mnm2m3 a m2 a m3 + F mnm2m . A yal n2 a m3 + ( ^o, m:i cil n2 J(15) 

Cmn = 2^6 i 1 + V") <5mn ' ^ mnw = 8^iv7^p' 5m+p+9 ' n (16) 
^ m " w = ^^^==S m , P+q+ n, F mnpq = 8 ^ §N ^__5 m+M+n (17) 

where repeated indices are always summed. Since our Hamiltonian is now in partially normal- 
ordered form, the Fock vacuum |0.) is manifestly an energy eigenstate with 

H.\0.) = Eo\0.). 

In the above computation we have retained all the commutator terms arising from reordering the 
a's and a' 's of the original Hamiltonian. Because our starting Hamiltonian was at most quartic in 
the fields the only such terms are c-numbers and quadratic. Furthermore, because of our treatment 
of p + , there are never reordering terms with only annihilation or only creation operators. However 
the reordering does alter the original color structure in the sense that we started with only a single 
trace and ended up with reordering terms involving two traces. Partial normal ordering of the 
general light-cone Hamiltonian is discussed in Sec. 5. 



3 Reduced momentum and reduced number operator 



In this section we begin our study of the large N matrix mechanics [16-18], [9], [19, 20] of light-cone 
field theory, drawing heavily on the methods and results of Ref.[9] 

We start with the simpler problems of the momentum operator P. = P + and the number 
operator N. 



P. = P f =Y, mTrala r . 



N. = N.^J2 Tra l 



N. 



a m \ 




- —m ( 


CL m 


VnJ t . s _ 






( a m \ 




= -( 


dm 


\VNj rs _ 







p.. 



N., 



N 



m 



N 



N 



P.\0.) = N.\0.) = {.0\P. 



(.0\N. = 



(18) 

(19) 

(20) 
(21) 



which are two examples of so-called trace class operators in the unreduced theory. In large iV 
matrix mechanics one then takes matrix elements of the commutation relations above, using only 
the ground state |0.) and the dominant time-independent adjoint eigenstates \rs,A) which saturate 
the traces of the theory at large N. The large iV completeness relation for these states reads (t = x + ) 



l. = \0.}(.0\ + J2\rs,A){rs,A\, 

rs.A 



d , > d . 



(22) 



because large N factorization guarantees that the singlet ground state will be sufficient to saturate 
the singlet channels of the traces. 
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Next one introduces reduced matrix elements^, for example 



(>#...£#) \st,B) 



|0.) 



N 



(.0|Tr 



Pl{h) Pn(tn) 



N 



N 



Pqp,rt (A\p 1 (t 1 )--- Pn {t n )\B) 
JV(0|pi(*i)---^ n (t n )|0) 



°sr,pq — 3rp&sq jy&sr&pqi 



a or a' 



(23) 
(24) 
(25) 



where U = xf and P is a Clebsch-Gordan coefficient. Matrix products (pi • • • p n )rs are called den- 
sities in the unreduced theory. These "Wigner-Eckart" relations then define the reduced operators 
a m and as well as the reduced ground state |0) and the reduced adjoint eigenstates \A). The 
reduced eigenstates satisfy the reduced completeness relation 



1 = |0><0|+X>><4, 

A 

at large N. 

The reduced forms of Eqs. (19) and (20) are 



d_ 

It 



|0> 



eft 



L4)=0 



(26) 



[P, a m ] 
N,a„ 



-ma r 



—a. 



P (V 

N,aL 



ma' 



(27) 
(28) 



Here P and N are the reduced momentum operator and the reduced number operator respectively, 
and the reduced operators a m and satisfy the Cuntz algebra 



O'mO'n — 



Smn, a ln a m = 1 ~ 1 0> <0| 

m 

a m \0) = (0\al = 



(29) 
(30) 



where |0) is the reduced ground state of the system. The Cuntz algebra is a simple example of a 
so-called free algebra, defined by products instead of commutators, and this free algebra is known 
to describe classical or Boltzmann statistics, another classical aspect of the large N limit. 

In fact the Cuntz algebra is a subalgebra of the so-called symmetric Cuntz algebra [9] which 
includes the following further operators and relations 



a m a{ = 5 mn , ^ aXfim = 1 - |0)(0| 

m 

m,al] = <W|0)(0|, [a m ,a n ] = [a^,a^] 
a m |0> = (0|a^ = 0, al\0)=al\0), (0|a m = <0|a, 



(31) 

(32) 
(33) 



Moreover, the tilde operators a m , aJ m are important in deriving the Cuntz algebra (29) from the 
unreduced theory: The reduction procedure (see Eq. (2.29) of Ref. [9]) gives only the relations in 
Eq. (32) as a consequence of the unreduced commutators (4), and a completeness argument (see 
Eq. (3.7) of Ref. [9]) on the general basis state is necessary to obtain the Cuntz algebra from (32) 
and the ground state conditions. Beyond this application, the tilde operators are not used in this 
paper. 

^Reduced matrix elements were introduced by K.Bardakci in Ref.[14]. 



5 



The reduction procedure does not give us directly the composite structure of reduced trace class 
operators such as P, N or the reduced Hamiltonian H. This is called the opacity phenomenon in 
Ref.[9]. But one can in principle find the composite structure of such operators by solving reduced 
commutation relations such as Eqs.(27) and (28). 

In fact, the reduced commutators (27) and (28) are so simple that standard methods [9] suffice 
to determine the form of the reduced operators P and N. To express these solutions succinctly, it 
is convenient to introduce a standard tool called word notation: 

(34) 
(35) 





— Q"m\ 


a m„ 


, w = mi ■ 


■m n 




= a \n n " 


"mi 


n 




M 


= n, 


{w} 







(36) 
i=i 

Here w is an arbitrary word composed of some number of ordered letters, which label the lattice 
sites in this case. In (36) we have also defined the length [w] of the word w and its weight {w}, 
which is the sum of its letters. 

Then the reduced momentum and reduced number operators can be written as 

P = P ] = a J I ma ln a m I a w 
w \ m / 

= ma^am + J2 a l\Yl ma L a m I «p H (37) 

m p \ m / 

TV = N ] = J2 a J {J2 a ln a m I a w 

w \ m / 

= Y a ™ a ™ +E4 E a m a ™ ) a p + • • • ( 38 ) 

m p \ m / 

P|0) = N\0) = (0\P = (0\N = (39) 

where J2 W m eans the sum over all words w, including the trivial word w = 0. To check that 
these operators satisfy the commutation relations (27) and (28), write them out, e.g. [P, a m ] = 
Pa m — a m P, and use the Cuntz algebra (29) to cancel all but one of the terms in the difference. 
Later we shall discuss more sophisticated ways to handle operators of this type. 



4 Reduced Hamiltonian 

The reduced momentum and number operators discussed above are very simple examples, but the 
opacity phenomenon has been particularly vexing in the case of models described by the equal-time 
quantized Hamiltonian 

H. = Tr ( -TT m 7r m + NV (^=) ) ■ (40) 



2 \^/N 
For such systems one knows the reduced equations of motion [9, 16] 



i[H, <p m ] = 7T m , TT m = i[H, 7T m ] = -V m ((f>) (41) 

p(t) = e iHt P (0)e- iHt , p = d> m or vr m (42) 
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where H is the reduced Hamiltonian, and one knows the interacting Cuntz algebras [9,19]. But 
beyond the case of matrix oscillators [9] and the general one-matrix model [17], the explicit form 
of the reduced Hamiltonian H of such systems is not known (an implicit construction of H is given 
in Subsec. 4.6 of Ref. [9]). 

For quantized matrix theories of the light-cone type, i.e. with no pure creation or pure annihila- 
tion terms in H., we shall see that we can in fact find the explicit form of the reduced Hamiltonian. 

As an example we begin with a general unreduced quartic Hamiltonian of the light-cone type 



H. — En 



A r , 



N 2 



N 



A r 



n' Vn , 



Cn 




(43) 



mnpq \ + 

p Q 



D 



mnpq i 



Cn 



D 



mnpq 



= c* n 

= D. 



N 

D, 



a p a q + 



Emnpq + . Fmnpq 



N 



mnp 



a 



nmqpi 



E 



mnpq 



nmpi 

= E, 



C, 



N 

~- F r , 



nmqpi 



mnpq - nmqp 

f>* D 



(44) 
(45) 

Here the quantities B,C, D, E, F are constants and the relations in (45) guarantee A\ nn = A nm 
and hence = H.. In what follows, we often refer to the form A mn in (43) as the kernel of the 
Hamiltonian. We emphasize that the general quartic form (43) includes the local quartic example 
in (12) and (13) as well as many light-cone models which are non-local in x~ (use Eq. (8)). The 
non-local models, including for example the number-conserving truncations of local models, can be 
considered as quantum-mechanical models in their own right. 

Following the standard procedure in large N matrix mechanics, we compute the time derivative 
of the annihilation operators 



TV 

(G m ) 



m) TS 



H. 



E 



{£m) rs = 



A 



rs 

mn^n 



r 



tu 



% Sir' q 13' 



rs-^mq 



N 



;Tr 



{Apq)i 



vt 



(46) 
(47) 

(48) 
(49) 



The term called T m comes from "opening" a trace in the double trace term of the Hamiltonian. 
To understand this term, we need the large N trace theorem [9, 17] 



N 



Trf 



n' Vn n 



<0|/(at,a)|0), V/ 



(50) 



which gives the leading term of the trace at large N in terms of the reduced average on the right. 
This implies in particular that 




iV' y/N , 







(51) 
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and so the term T m fails to contribute to the large TV equation of motion. 

The term called £ m collects the contribution of all creation operators except the first in the 
a^Aa term of the Hamiltonian. By direct computation, we find the explicit form of this term 



i (£ m ) r 



1 

N 

+Q 



a 



pqrn 



pqmx 
~\~-Dpqmx 

Each of these terms has the form 




pqxm 



pqxm 



4 




a q 


y/N 






4 


a x \ / 


' a q 


y/N 







rt 



vt 



(52) 



called a "twisted density" in Ref.[9]. Twisted densities can be reduced directly in terms of the 
tilde operators in Eq. (33), but we prefer to express the result first in terms of untwisted densities 
(ABCD) rs . Thus we commute the factors into the form F rt E ts — > (FE) rs , keeping all correction 
terms. 

The result of this algebra is 




nVn, 



Cppm^rs 



a q dp 



y/N 


y/N 


y/N t 


a q 


4 


a x 


y/N 


y/N 


y/N \ 


a x 


a q 


4 


y/N 


y/N 


y/N 


a q 


4 


a x 



N y/N y/N , 



ts 




1 at 

C, A __Tr 1 

y- / pqqm u rs jy - LL 



D 



ppmx 



J ppqm 



(53) 



where each term in parentheses corresponds to one of the twisted densities in Eq.(52). The two 
trace terms in (53) can again be neglected by the large N trace theorem (51). 

Now all surviving terms in the large N equation of motion are densities, so we may use reduced 
matrix elements to obtain the reduced equation of motion 



ia r 



= [H, a m ] = A mn (a\a)a n + Extra (54) 

Extra = {^CpqmOiqc/p Cpprnj ~\~ {/-^pqmxOj x (lq0^p CpprnxQ^/j ~\~ {CpqxmQ-qt/p^x Cppxm4) 

~l~ ^Dp qmx a, x <i q (i}p DppmxCi/j + (^Ep qxm a q a^p(i x — Eppq m ciq*j (55) 



where A mn (a) , a) is the reduced form of the kernel, given explicitly below. The parentheses here 
correspond to the parentheses of Eq.(53). But now we notice that each of the extra terms in 
parentheses is separately zero, e.g. 



EpqxmO'qo)pCL x — Epp Xm d x 



(56) 
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by the Cuntz algebra. 

The results of our computation are the simple reduced equations of motion 

a m = i[H,a m ] = -iA mn (a\a)a n 

a L = i[H,a\ n ] = ia\ l A nm {a ] ,a) (57) 

where H is the reduced Hamiltonian. 

For such simple reduced equations of motion, standard methods again suffice to write down the 
reduced Hamiltonian of the large N light-cone system 

H - E = J2 a wHJ2 a ln A rnn(a t ,a)a n )a w (58) 

w m,n 

= ^a\ n A mn (a\a)a n + ^al(^al n A mn (a\a)a n )a p ^ (59) 

m,n p m,n 

A mn ((l , 0") — C mn + Cmnp&p ~\~ E , mn p(lp + CmnpqOpO^q 

~\~E) mn pqa)paq + E mn pqapQ^q + F rnn pqapCLq (60) 

(H - E )\0) = (0\(H - E ) = (61) 

where the form A mn (a\a) is just the kernel we started with in the unreduced Hamiltonian (43). 
The reduced Hamiltonian is hermitian because A] nn = A nm . The Cuntz algebra also tells us that 
the E term of A mn can be further simplified 

EmnpqttpO^q = E mn pp (62) 

so that the kernel A mn can be considered as completely normal-ordered. 

In the following section we will argue that the results (57) and (58) are in fact the forms obtained 
by the large N reduction of any light-cone system, where the reduced kernel A mn is always defined 
by the single trace term of the Hamiltonian. As a check on this conclusion, we note here that the 
reduced equations of motion (57) are consistent for all A mn with the ground state conditions (30), 
for example 

a m |0) = -iA mn (a f ,a)a n \0) = 0. (63) 
Moreover the reduced equations of motion are consistent with the Cuntz algebra (29) 



d 

-r I W 

at 



n) = |(E t 4am)=0 (64) 



for all A mn . Note that the Cuntz algebra itself is necessary 

(a m a+) = i(-A m p(a\a)a. p a\ l + a m a^ p „(a t , a)) = (65) 
to check the first of these relations from (57). 



5 The General Reduced Light- Cone Hamiltonian 

Now we extend this discussion to a general light-cone Hamiltonian. First consider a Hamiltonian, 
not necessarily of the light-cone type, required to be a color singlet and behaving at large iV no 
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worse than 0(N 2 ). To implement these restrictions we introduce a Fock space, which can always 
be done, but in the general non-light-cone situation the Fock vacuum is not an energy eigenstate. 
To count powers of N, we first recall the well known estimate 



(.0|TrF^,^|0.) = O(AT) (66) 



with F an arbitrary polynomial of its arguments. This estimate is a simple consequence of the 
matrix nature of the operators and Wick's theorem. The fact that the vacuum expectation value 
of a single trace gives the maximal power of N is embodied in the large N factorization theorem 

(.0|TriqTrF 2 |0.) ~ (.0|TrFi|0.)(.0|TrF 2 |0.) = 0(N 2 ). (67) 

Thus a generic color singlet Hamiltonian scaling as N 2 will have the form 

H = N 2 e + NTrhn + Trh 21 Tih 22 + ■■■ + j^Tvh kl TTh k2 ■ ■ ■ Trh kk + ■■■ (68) 

where eo is a c-number and each h is a matrix function of a/y/N and a) /y/N. Without loss of 
generality, we may put the operators in each trace in a partially normal-ordered form as explained 
in Section 2. This is done by always moving an annihilation operator at the extreme left of a trace to 
the extreme right or a creation operator on the extreme right to the extreme left. The commutator 
terms have a creation operator and an annihilation operator removed and an additional trace times 
1/N. They therefore preserve the generic structure of (68), simply redefining the h k i- This process 
is continued until every trace has the desired form. Note that this reordering procedure preserves 
the cyclic ordering within each trace, but modifies the h^s in terms with a smaller number of a's 
and o''s. Moreover, we can assume that each h k i has zero vacuum expectation value: for if it does 
not then the VEV has the right N dependence so that its contribution can be cancelled by simply 
modifying one of the other generic terms. 

Now we describe the further restrictions enjoyed by light-cone Hamiltonians. In this case every 
operator term in the Hamiltonian has at least one annihilation operator and at least one creation 
operator. In our generic Hamiltonian (68) we can therefore order things so that in every term 
there is a creation operator on the extreme left and an annihilation operator on the extreme right. 
With this ordering prescription |0.) is manifestly an eigenstate of H with eigenvalue Eq = N 2 eo. 
Although |0.) is not necessarily the lowest energy eigenstate, we shall assume in what follows that 
it is. In more detail, our ordering prescription shows that the /i's have the form 

t f 

hi i = — Ad ; Hua — — ^=Bu i ; huu = CVt. — ■== ; hhi = — 1= J3i.j + CV/ — ^= (69) 
TV Jn y/N \/N vn 



for k > 1 and 1 < I < k. Here we suppress all indices on the operators. Note that A in h\\ is 
the kernel of the Hamiltonian defined in Sec. 4, and we remark that Eq.(68) with the restrictions 
of (69) generalizes the discussion of the quartic light-cone field theories of Sections 2 and 4 to an 
arbitrary Hamiltonian of the light-cone type. 

The first step in inferring the dynamics of the reduced formulation is to work out the equations 
of motion for the unreduced a rs 's. The Heisenberg equation for an annihilation operator is 



at 



y/N \ \y/N' y/N J y/N 



y/N ,A<y y/N 7 y/N^ 



a + F (70) 



T = U WTr() + -^( ) rs Tr()Tr() + ..., (71) 
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where we have suppressed the arguments of the various densities and traces in T . 

The equation of motion for the reduced annihilation operator should coincide with the unreduced 
equations of motion at N — > oo for matrix elements between the singlet and dominant adjoint 
eigenstates as in Eqs. (23) and (24). Since the dependence on the open color indices is completely 
determined by group theory, it is sufficient to test the coincidence of equations of motion in the 
reduced and unreduced formulations for these matrix elements with all the open indices in Eqs. (23) 
and (24) saturated. This is equivalent to testing their coincidence when the unreduced equation of 
motion is placed within single trace VEV's of the form 

(.0\Ti(Df(t)E)\0.), VD,E (72) 

where f(t) = is the unreduced equation of motion and where D and E are arbitrary matrix 
functions of the unreduced a's and a^'s. Any simplifications at large N in f(t) valid within all such 
single trace VEV's can be exploited in deducing the reduced equations of motion. 

We next note the simplifications that occur when the right side of (71) is inserted in such a VEV. 
Although every term is nominally of the same order, the multi-trace terms T are all negligible at 
infinite N because, by large N factorization, each trace can be replaced by its vacuum expectation 
value, Tr(hki) — > (.0|Tr(/i^)|0.). And the VEV of each hya vanishes by our ordering convention. 
Thus all the terms of T in (71) can be immediately dropped. 

The commutator in the second term on the right side of (71) is a sum of terms of the form 



Tra 1 " 



""=. L(4t>4t) a~^{cJB) us {Ca) ru . (73) 



y/N' WN' y/N 



each of which is a twisted density. These are the terms which were most tedious in the example 
of Sec. 4. We will show that, when inserted in the VEV of a single trace, such twisted density 
terms are always suppressed relative to nominal order by two powers of 1/N. Then the right side 
of the equation of motion can be replaced by only the first term, and the reduction procedure is 
straightforward. 

First we confirm that the first term can contribute at nominal order. So consider 

{.0\Dtr (A-^) E st \0.), VD,E. (74) 

The explicit a to the right of A will contract against an operator in E. If E has a creation operator 
on the left, i.e. if E st = a\ u E' ut / 'y/N \ the contraction with this operator gives a factor 5 SS /N = 1 
so the resulting term will give 

(.0\D tr A rs E' st \0.) = (.0\Tr(DAE')\0.) = O(N). (75) 

Of course, whether this leading term is actually present depends on the detailed structure of DAE, 
e.g. (.OlTraValO.) = even though it is of nominal order N. 

Next we give the corresponding evaluation for the twisted density terms on the right side of 
Eq. (71): 

(■0\D tr (^=B) us (C-^=) ru E st \0.). (76) 
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In this case the a on the right will contract only against operators in E and the a) on the left 
will contract only against operators in D: 

(■0\Dtr(—7=B) us (C—7=) ru E st \0.) = {■0\D tr (—7=B) us — (C) rv F su G v t\0.) 
V-/V y/Jy y/Jy iv 



= (.0\HtaK ur ^B as ^(C) rv F su G vt \0.). (77) 

Following the color indices, it is easy to see that this involves a single trace, albeit in a "twisted 
form" . But this single trace is multiplied by the further factor l/N 2 , so the contribution is nominally 
of 0(1/N). This is suppressed by two powers of N compared to the leading contribution (75). We 
conclude that 

<.0|Tr(D (id(t)) E)\0.) = (.0|Tr(D (Aa) E)\Q.) (78) 

for N — > oo for arbitrary matrix operators D, E. In the language of the previous section we have 
shown that 

l^B) us (Ca) ru = (79) 

^177.^177 (80) 



when sandwiched between the ground state and the dominant adjoint eigenstates at large N. Since 
this simple large N equation of motion involves only densities, it is easily translated into the reduced 
equation of motion given in (57) 



la 



= A(a),a)a (81) 



and one finds that the reduced Hamiltonian of the system is that given in (58), now for arbitrary 
reduced kernel A. 

Finally we must show that this result does not depend on our choice of partial normal order- 
ing. Indeed, the different ordering decisions lead to apparently different expressions for the kernel 
A (^a) , aj . These differences occur in the terms in A that are not completely normal ordered. How- 
ever, as noted in Sec. 4, use of the Cuntz algebra shows that these terms immediately collapse to 
completely normal-ordered ones and the discrepancies disappear. 

This last point can also be seen in the unreduced formulation. Consider the contribution of a 
term in A of the form Ba(x)a(y)^C/N to the right side of Eq. (78) 

— {.0\Tr(DBa(x)at(y)CaE)\0.) = (.0\Ti(DBCaE)\0.)c(x,y) 

+ X) ^<-°l[ £>S ]«' :i> [ at (2') C "][ Cf " ajE7 ]«-l -> + ^<-0|[J3S]„TV[at(j / )c a .E?'][£?"] w |0.> 
~ (.0|Tr(D£Ca£)|0.)^ 
+ J2 {.0\ r Tr[DBE"]\0.){.0\Tr[a t {y)CaE , ]\0.) = (.0\Tr (DBCaE)\0.)c(x, y). (82) 

Here we have used x, y to collectively denote all the labels carried by a, a^. Then we defined c(x, y) 
by [a rs (x), a\u(y)] = Sst8ruc(x,y). Thus at infinite N, wherever we find the combination a(x)a f (y) 
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in a term in A we can replace it by the c-number c(x, y). This is a mirror in the unreduced theory 
of the action of the Cuntz algebra in the reduced theory: After this replacement has been made 
wherever possible, every term in the simplified A will be in the completely normal-ordered form 
a) n a m . This simplified form of A, which we call A^, is then unique. In summary, for the purposes 
of the large N limit we may replace the unreduced Hamiltonian H. by 

H. = Hqo = Travel (83) 
N 

where Aoo is the completely normal-ordered kernel defined above. 



6 The Algebra of All Normal- Ordered Operators 

Having completed the reduction of the general light-cone theory, we discuss in this and the follow- 
ing sections some simple algebraic aspects of the reduced Hamiltonian systems (58) for arbitrary 
normal-ordered kernel A. In particular, the infinite dimensional algebra developed from the Cuntz 
algebra in this section is a simpler version of the infinite dimensional algebra developed from the 
interacting Cuntz algebras in Ref. [19]. 

We begin by expressing the Cuntz algebra as an infinite dimensional free algebra [19] 

a w a w ^ = 5 WtW i + y^(5 WtUW >a u + 8 w i tUW aJ) (84) 

a w a w > = a ww >, aja w ^ = a w ,J (85) 

x f 1 when w = w' toa \ 

S w , w > = < . (86) 

I U otherwise 

where w, w' and u are arbitrary words. 

Next we introduce the general normal-ordered product of Cuntz operators 

E w - W > = aja w >, E W ; W > = E W;w, Vw,w' (87) 
and we use (84,85) to obtain the infinite dimensional algebra of all normal-ordered Cuntz operators: 

E U \u' E v;v' ^u' ,v E u;v' ^ ] i^u 1 ,wv E u;wv' $v,wu' E wu;v'^) 

= ^ ' L uu i - vv i ; X x' E x;x' ($8) 
x,x' 

Luu';vv';xx' = $u' ,v^u,x^v' ,x' ^ ' i&u' ,wv$u,x$wv' ,x' ^v,wu'^wu,x^v' ,x') ($9) 

\_Eu;u' i E v;v'1 ^ ' E uu i- vv i- xx i E x - X i i E uu' ;vv' ;xx' I J uu';vv';xx' ^vv' ;uu' ;xx' ■ (90) 

x,x' 

Here x and x' are also arbitrary words. In spite of appearances, the right sides of the product 
relations (84) or (88) yield exactly one term, for example 

E U . U ,E V . V > = E u . wx , when u' = wv,w^0. (91) 

We also note that 

[Eo; W , Eyj/.fi] = 5 WyW ' — E w i. w + ($w,w"w'Eq- w » + 6 w i, w h w E w ».q) (92) 
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where 5 W w i is a central term proportional to Eq-q = 1 which commutes with all the generators 
of the infinite dimensional algebra. Subalgebras of the algebra of normal-ordered operators are 
considered in the following section. 

Our reduced number and momentum operators N and P are easily written in terms of these 
operators 



N 



^ ' Q"w I ^ ' a rn®"ni ) a w — ^ ] E mw - jT 
w \ m ) mw 



p = E^ f E 



m / mw 



(93) 
(94) 



Moreover for a general normal-ordered kernel A mn we may also express the general reduced Hamil- 
tonian = H as 



Ct 



CI A^ - A 

^nmvui -"-ran — ^nm 



H - E = ^aJ(^a\ n A mn (a\a)a n )a w = ^ C mnuv E,, 



umw;vnw 



(95) 
(96) 
(97) 



w m,n 



w,u,v,m,n 



where u,v and w are arbitrary words and C mnuv are constants. 

Then we may reconsider various commutators above as consequences of the infinite dimensional 
algebra (90). For example it is not difficult to check that 



N,a m 
N,al 
[P, a m ] 



P of 

1 i u m 



^ ' [Enw;nwi ^0;m] — -^0;m — 
nw 

^ ] [Enw;nwi E m -o\ = E m -Q = o) m 

nw 

^ ' 72 [E nw - nw , £/0;m] = TTlCt m 
nw 

E n [Enw;nw, E m -o] = TXIO^ 



(98) 
(99) 



[H, a m ] = 



A m nQ-ni 



H,al 



dn-^-rim- 



From (98,99) we find also that the reduced P commutes with all zero momentum operators, e.g. 

P, a m a n o)p(i m ^ n j t .p = [P, Ep nm - >m + n +p\ = 0. 
The general form of this statement is 



[P,E W . W >] = {{w} - {w'})E w . w > 

[P, E w . w i] = when {w} = {w'} 



(100) 
(101) 



where {w} is the weight of w (see Eq. (36)). It also follows that the reduced momentum operator 
is conserved 

[P,H}=0 
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when the kernel A mn and hence H is constructed from zero momentum operators. Similarly the 
reduced number operator commutes with all zero number operators 



N E i 



(M - [w'])E w . w , 
when [w\ = [w'], 



(102) 
(103) 



where [w] is the length of w (see Eq. (36)). Then the reduced number operator always commutes 
with the reduced momentum operator 



N,P 







(104) 



and the reduced number operator is conserved [H,N] = when A mn and hence H is composed 
of zero number operators. We shall return to number-conserving reduced Hamiltonians in the 
following section. 

We also consider the action of these operators on a general basis state 



H = aj\0) = a) mn •••aj ni |0>, Hm) = l 
and for this discussion the following identities are useful 

E U ;v\w) = ^2S WjW > v \w'u) 
w' 

^ 5 w ,w'mw" = [w] 
w',m,w" 

mS w,w'mw" = {w}. 

w' ,m,w" 

Here [w] and {w} are the length and weight of w, and we find that 

N\w) = [w]\w), P\w) = {w}\w} 
as expected. The action of the reduced Hamiltonian (97) on the general basis state is 
(H - E )\w) = ^2 C mnuv 5 WtW > vnw// \w'umw"). 

m,n,u,v,w',w" 

Similarly, one may evaluate the general matrix element of H — Eq 



(w'\ = (Q\a w i, {w'\w) = 5 w i. 

■"mnuv$w,w"'vnw" u w' ,w"'umw" 



(w'\(H — Eq)\w) — C m nuv5w.w"'imw"Sw'.-'"' .,,,,„■" 



m,n,u,v,w",w" 



and so on. 

As a simple illustration, we mention the one-matrix cubic interaction: 



H = J2a^ n Aa n , A = cu + A(o + a f ), \n) = a tn |0) 



n=l 



(H - E )\n) = nu Q \n) + \(n\n+l) + (n-l)\n-l)), n > 0. 



(105) 

(106) 
(107) 
(108) 

(109) 
(110) 



'111) 



(112) 
(113) 



We solved this simple recursion relation numerically, but this solution is superceded by a recent 
analytic solution given in Ref. [22]. 
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We finally note that, in large N matrix mechanics, the large N master field [23] is identified 
[16] as the set of reduced Heisenberg operators in an energy eigenbasis. For example, we have 

(H-E )\a) = E a \a), a = (0,A) (114) 

a m (t) = e mt a m {Q)e~ mt (115) 

a m (t) a p = e^o m (0) Q/3 , uj a p = E a -Ep (116) 

a{t) a p = iaj al3 a(t) a p = -iA(t) ai a(t) l/3 (117) 

so that the time dependence of the master field is controlled by the energy differences w Q/ g among the 
singlet ground state |0) and the dominant adjoint eigenstates \A). These are the time-independent 
states which appear in the reduced completeness relation (26). We return to eigenvalue problems 
in Sections 7 and 9. 

7 Subalgebras and a Class of Integrable Models at Large N 

We begin this section with a partial list of subalgebras of the infinite dimensional commutator 
algebra (90) of normal-ordered operators. 

1. The product algebra (88), and hence the commutator algebra (90), of all zero momentum 
operators 

E w ,w', {w} = {w'} 

are closed subalgebras. Proof: Consider the product of two such operators, noting that the 
Kronecker deltas on the right of (88) always pick out a unique operator - which is then a 
zero momentum operator. (Or note that P commutes with a product of two zero momentum 
operators.) 

2. The product and commutator algebras of all zero number operators 

E w>w >, [w] = [w'] 

are closed subalgebras. The proof follows as above for zero momentum operators. In this 
case we list some examples 

[Prn;ri7 Pp;q\ — &npPm;q &m,qPp;n (H^) 
\Emn;pqi Pl;r] = dqlPmn;pr ~ ^nrPml;pq 
{Pmn;pqi Plr;st] = 3qr&plPmn;st ^nt^msPlr;pq (1^0) 

with m, n,p, q = 1 . . . K. In particular we recognize (118) as the algebra of the unitary group 
U(K), where K is the number of lattice sites. 

3. Cartan subalgebra and subalgebras thereof. These can be useful in finding integrable systems. 

It is not difficult to find one subalgebra of the Cartan subalgebra of (90). Consider the hermitian 
"single word" operators 

Pw;w — "^"ui! E W ' W ^ — E w - tWl \/lV 

which are a subset of the zero momentum and zero number operators. These satisfy the product 
relation 

Pu;u-Ev;v = &u,vPu;v + ^ ' {^u,wvPu;u + $v,wuPv;v) (121) 



16 



and hence each single word operator E u - U is a projection operator 

(E U]U ) 2 = E u . u . (122) 
Moreover the relation among the projection operators 

{E u;u when u = wv 
E v . v when v = wu (123) 
otherwise 

is an equivalent form of Eq. (121). Finally, we find that all the single word operators commute 

[E U ; U ,E V . V ] = 0, Vu,v (124) 

because the right side of (121) is symmetric under u v. 

This leads us to consider the following large class of momentum and number-conserving reduced 
Hamiltonians 

^mn = 0~mn ^""^ X W E W - W (125) 

w 

H — Eq = y ' X w E wjnw i. wmw i = ^ ' ^w a wmw ' a wmw' (126) 



wmw' wmw' 



where {X w } is a set of arbitrary word-dependent couplings. It follows from (124) that each of these 
models possesses an infinite number of conserved quantities 



d d i | 

Jt w ' w = Jt\ aw 



ajj = 0, Vu> (127) 



which leads us to expect that each of the Hamiltonians (126) is in fact an integrable model. 

As further evidence that these models are integrable, we have checked that their eigenvalue 
problems are diagonal on every basis state \w): 

(H - E )\w) = E(w)\w), E(w)=^2\ w f ^2 S w,w"'w'mw" ■ (128) 

w' w"'mw" 

We can be more explicit by choosing a simple word dependence A^ = A ([to]) for the couplings and 
using the identities 

]T 5 WjW , w „ = [w] + 1 (129) 

w',w" 

S w,w'mw" = N (130) 



w ,m,w" 



sr a _ i H - i H > 2 

<Jw,w l mnw" — 1 Q [to] = 1 \ 101 ) 



w'm,n,w" 



E 

letters of u at fixed [u] 



w 



M + i 



for [w] > [u] 
for [w] = 0, 1, . 



u 



(132) 
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which generalize Eq. (107). In this case we find the eigenvalues 



M-i 

E(w) = ]T A(M)(H - [«]). (133) 

[u]=0 

Note that each term A ([it]) shows a shifted oscillator spectrum [w] — [u]. 

This class of integrable models is very special. As number-conserving interactions, they corre- 
spond to unreduced Hamiltonians which are non-local in x~, but moreover these particular models 
provide only a subset of the number-conserving terms of physical light-cone models 



1 f f 

-. ■ r a m a n a p a m+n-p 

m n y / mnp{m + n - p) 



5 j (^m^n^n^m — 5 1 E nm - nm . (134) 
p =n m,n m,n 



We emphasize however that other integrable models would follow should the Cartan subalgebra of 
the infinite dimensional algebra (90) prove to be larger than the subalgebra (124). 



8 Hidden Conserved Quantities at Large N 

In Ref.[9] it was pointed out that reduced formulations show new hidden conserved quantities at 
large ./V and, moreover, that these can be pulled back to find hidden conserved densities at large 
N in the unreduced theories. These features are shared by our reduced light-cone models. 

In the first place, the Cuntz algebra itself provides us with a set of new conserved quantities 

[a m a\^ = a tn a mJ = 0, m,n = l...K (135) 

for all models in our class. According to the reduction procedure, these conserved quantities can 
be pulled back to find new hidden conserved densities in the unreduced theory 



dt\ N N dt\„ N N 



rs 



The notation here means that these densities will be conserved at large N, at least when sandwiched 
between the dominant unreduced eigenstates |0.) and \rs, A) which saturate the traces of the theory. 

Furthermore, every time we succeed in finding the explicit form of a reduced trace class operator, 
we get another hidden conserved density in the unreduced theory at large N. In our case we find 
the new densities 



dt r " N 



= o, p rs ^Y, N ' H ~ 1 \ a w t \T, ma l a ^) a ^) ( 137 ) 

N w \ \ m / J rs 

jH rs = 0, i? rs e ^ JV-N- 1 ^ aUmn ^] a n a w j (138) 



which correspond to the explicit forms (37) and (97) of the reduced P and H. Here the words are 
constructed as matrix products from the matrix creation and annihilation operators 

(fio) rs = (Omi ' ' 1 fm„) rs ) \ a w) = ( fl L„"' a Li) 
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Following Ref.[9], we remark that the first terms of the traces of these densities are proportional to 
the dominant parts of the original unreduced trace class operators P. and H. 



TrP 



Jj R + ]^2 Tr 2 (4 (ll ma ln a rr)j <hj + ' ' ' ( 139 ) 

^(E«L^„(^,^) «„)+.-. (140) 



M = iV 

but there are an infinite number of other terms from the "dressing" of these densities. 

Finally, for each of the integrable models of Sec. 7, we find that the entire infinite set (127) of 
reduced constants of the motion can be pulled back 

±( N -^aJa w ) = 0, Vw (141) 
at \ / rs TV 

to appear as hidden conserved densities in the unreduced theory at large N. 



9 Extensions 

9.1 The dominant states in the unreduced formulation 

Although we have reformulated the large N eigenvalue problem in a simpler free-algebraic form, 
it is also useful to identify the basis of the dominant adjoint eigenstates \rs, A) in the unreduced 
theory. To do this we return to the discussion of Sec. 5. 

In particular, we can insert the unreduced Hamiltonian itself deep inside the VEV of a trace. 
By large N factorization the leading behavior would be 

(.0\Tt{DH.E)\0.) ~ (.0\Tt(DE)\0.)(.0\H.\0.) = 0{N A ) (142) 

since the Hamiltonian is itself a trace. The interesting physics, however, is contained in the con- 
nected part of the VEV which is of order N 2 . We can isolate the connected part by performing the 
contractions of an explicit annihilation operator in H. with operators in E or those of an explicit 
creation operator in H. with operators in D. Again all the terms in H. with two or more traces 
(see Eq. (68)) are negligible and we obtain 

(.0\Tt{DH.E)\0.) - (.0\Tt{DE)\0.}(.0\H.\0.) ~ (.0|Tr (zTTr (oUa) E) |0.) c (143) 

where the subscript c indicates that at least one of the contractions noted above is taken. We are 
thus led to consider the following action of H. on the special basis states 

|0.) (144) 

rs 

as was done for singlet states in [11], where it is shown that nearest neighbor contractions dominate 
at infinite N. The extra \fN in (144) normalizes these basis states. This is an alternative path to 
the Cuntz form of the reduced Hamiltonian at large N. It is also noteworthy that the content of the 
Cuntz algebra in the reduced space is essentially summarized by this nearest neighbor interaction 
rule. The special basis states in (144) are in 1-1 correspondence with the reduced basis states (105), 
and the dominant adjoint eigenstates \rs, A) are linear combinations of the basis states in (144). 
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9.2 Energy eigenvalue problem in the singlet sector 

Although we do not have a free-algebraic formulation of the singlet sector of our matrix models, 
similar simplifications have been found in the corresponding unreduced large N singlet eigenvalue 
problems [11]. For example, we find in the case of the integrable models of Sec. 7: 



H.(u) = A^-MTr (a u f a u ) (145) 

\w.) = 7V-M/ 2 Tr ( a j) |0.), (.w\w.) = l (146) 

H.(u)\w.) = E u (w)\w.), E{u)= J2 ( 147 ) 

w^vuv' 



where w and u are arbitrary words. Multiple trace terms which arise in this computation are 
negligible at large iV [11]. The notation used in the eigenvalue E u (w) was defined in Ref.[9]: Sum 
over all words v and v' such that vuv' is equivalent to w up to cyclic permutation of their letters. 

9.3 Light-cone Hamiltonian for Non-Abelian Gauge Theory 

When the methods of this paper are applied to Yang-Mills theories, the unreduced Hamiltonian is 
much more complex. In light-cone gauge, A_ = the physical potentials are 

A(x,x") = ^ 00 -^=(a l (p + ,x)e—^ + +aJ(p + ,x)e + »^ + ) (148) 

E^(«M(x)e-^ m5 + ai/x)e + ^^) (149) 
l=i v4tt/ 

where we have discretized p + in the last line. We see that the operators carry a polarization index 
i = 1, 2 as well as a transverse coordinate x, which can also be discretized. If one formally eliminates 
the longitudinal "Coulomb" potential A + , the unreduced Hamiltonian takes the form 

r i i 

H. = dyidx' -TrdiAjdiAj + igTxdiAj[Ai, Aj] + igTrdiAi— [Aj, d^Aj] 

J £ (J — 

- 9 ^Tr AiAj^Aj] - ^Tr (J-fa, d-A^ (150) 

r l i 

= / dyidx' -TrdiAjdiAj + 2igTrAj — [d-A u djAi] 

- 9 ^Tr AiAj^Aj] - ^Tr (j-[ Al , d-A^ . (151) 

Here the precise definition of the inverse derivatives dZ 1 has been left unspecified. 

The terms arising from partial normal ordering will all be of the form of a quadratically divergent 
gluon mass term. There must be counter terms that cancel not only this mass, but also contributions 
to the mass that arise in higher orders of perturbation theory. We must also allow for counter terms 
which remove potential Lorentz invariance violations induced by ultraviolet divergences in this 
noncovariant gauge. Such counter terms need not be local in x" . As a measure of the complexity of 
this Hamiltonian, we remark that in 4 space-time dimensions there are 12 cubic terms corresponding 
to different spin transitions, as shown for example in [24]. There are clearly many technical issues 
that must be resolved before taking this Hamiltonian completely seriously. 
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Even so, it is clear that once H. is completely specified, the results of this paper show that the 
reduced Hamiltonian has the Cuntz-algebraic form discussed above. Burying all of the unresolved 
issues in the kernel A, we obtain a reduced Cuntz-algebraic system of the form 



H — Eq 



E^ f ( E 



<X W 

E a lm( X K-( X ) = 1 -|0)(0| 



(152) 



«i,m(x)a] n (y) 




(153) 



a iim (x)|0) = (0|at m (y)=0, 



(H-Eq)\0) = (0\{H-E o ) = (154) 



where |0) is the reduced ground state. Moreover, the question of confinement can be addressed by 
looking for reduced eigenstates 



where W(x, y) is the set of all words with the transverse coordinates of the first and last letters 
constrained to be x and y respectively. For example, if we could find the lowest energy eigenstate 
subject to the constraint y — x = R, then linear growth of this energy with R = |R| would signal 
confinement [25]. 
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